Plasmon excitations in homogeneous neutron star matter. 
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ABSTRACT 

O ; 

We study the possible collective plasma modes which can affect neutron-star ther- 
modynamics and different elementary processes in the baryonic density range be- 
>- ■ tween nuclear saturation (p ) and 3p - In this region, the expected constituents 

(N : 

. of neutron-star matter are mainly neutrons, protons, electrons and muons (npe/i 

m ■ 

matter), under the constraint of beta equilibrium. The elementary plasma ex- 
citations of the pen three-fluid medium are studied in the RPA framework. We 
emphasize the relevance of the Coulomb interaction among the three species, in 



particular the interplay of the electron and muon screening in suppressing the 
possible proton plasma mode, which is converted into a sound-like mode. The 
Coulomb interaction alone is able to produce a variety of excitation branches and 
the full spectral function shows a rich structure at different energy. The genuine 
plasmon mode is pushed at high energy and it contains mainly an electron com- 
ponent with a substantial muon component, which increases with density. The 
plasmon is undamped for not too large momentum and is expected to be hardly 
affected by the nuclear interaction. All the other branches, which fall below the 
plasmon, are damped or over-damped. 

PACS : 21.65.+f , 24.10. Cn , 26.60. +c , 03.75.Ss 



I. INTRODUCTION 



The evolution of neutron stars is in general determined by the microscopic 
processes that can occur in their interior and in their crust. In particular, the 
cooling process by neutrino emission is strongly affected by the detailed structure 
of the matter and the excitations it can sustain. Mean free path and emissivity 
of neutrinos are directly related to the correlations present in the neutron star 
matter, which in the standard models is composed of neutrons, protons, electron 
and muons [1]. The effect of collective excitations on neutrino emissivity have 
been extensively studied in the literature, often with controversial results [2-8]. 
The proton plasmon excitations seem to play a relevant role in quenching the 
contribution of the pair-breaking mode to neutrino emissivity [5] by charged cur- 
rent, since it obscures the possible effect of the Goldstone mode, which should 
be present in a superfluid due to gauge invariance. However at low momentum 
the electron and muon screening is quite effective and we will show that indeed 
the proton plasmon does not exist, and therefore its role in neutrino emissivity is 
quite delicate and must be examined with care. For the same reasons, at low mo- 
menta the proton oscillations are accompanied necessarily by electron and muon 
oscillations, and this can be of great relevance on the total neutrino emission [7] . 
The real plasmon mode in neutron star matter is at higher energy and we will 
show that it is mainly an electron excitation, with very small contribution of the 
proton motion but a no n- negligible muon component. 

All that can be best studied considering the proton, electron and muon spectral 
functions, that can be used also to determine the damping of each mode. We will 
present extensive results on the spectral functions of the electron-muon-proton 
plasma under neutron star physical conditions. Full spectrum and width of the 
excitations will be presented in detail. 

Since neutron and proton are coupled by nuclear forces, it is not clear to what 
extent the presence of the plasmon mode can affect the whole set of collective ex- 
citations or the excitation spectrum at a qualitative level, even for non-superfluid 
matter. Some considerations and preliminary results will be presented on this 
point, which will stress that the genuine plasma excitations are in general ap- 
proximately decoupled from the rest of the spectrum, but in particular cases 
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they can be strongly damped. 



II. THE ELECTRON-MUON-PROTON PLASMA 

We consider protons, muons and electrons, in neutron star matter conditions, 
interacting by Coulomb coupling, disregarding the nuclear interaction of protons 
with neutrons. This idealized case, which of course does not correspond to actual 
physical conditions, will permit to clarify two points, also relevant for the general 
physical situation. On one hand it will make more transparent the role of electron 
and muon screening on the proton dynamics, on the other hand it will give clear 
indications on the strengths of the electrons, muons and protons motion at the 
different excitation energies. The role of nuclear forces will be briefly discussed 
at the end of the Section. 

As already mentioned in the Introduction, we are considering normal nuclear 
matter, leaving the case of superfluid matter to a future work. The range of total 
baryon density pb will be fixed between the saturation density p = 0.16 fm -3 
and 3po, where indeed neutron star matter is expected to be homogeneous and 
not affected by possible "exotic" components like hyperons. Of course the proton 
fraction Y p = p v j pb has to be fixed in order to choose the proton density, which 
is actually the only relevant parameter for the system under study. In order to 
specify Y p according to neutron-star structure, it is determined by beta equilib- 
rium, which depends on the nuclear matter EOS. This will be taken from the 
microscopic calculations of Ref. [9, 10], based on the many-body theory within 
the Bethe-Brueckner-Goldstone (BBG) scheme. For simplicity we will work at 
zero temperature, but the results can be easily extended to finite temperature. 
Electroneutrality imposes p p = p e + Pfi, and the proportion between electrons and 
muons is fixed by chemical equilibrium. The resulting proton and muon fractions 
as a function of baryon density are reported in Fig. 1. On this figure, the micro- 
scopic calculations are compared with results from three modern Skyrme forces: 
SLy230a [11], NRAPR [12] and LNS [13], in order to give the order of magnitude 
of the dependence on the EOS. The particle fractions used for the three baryonic 
densities under study are listed in Tab. I. 

For a one component plasma, i.e. a charged Fermi gas on a rigid uniform 
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background (jelly model), the Random Phase Approximation (RPA) has been 
studied in detail, both in the non-relativistic [14] and relativistic [15] cases. The 
excitation spectrum, i.e. the line in the energy-momentum plain where the real 
part of the polarization function has poles, is known to have a "thumb like" shape. 
This is illustrated in Fig. 2 (full line), where the case of a proton gas at a density 
p = po is considered. 

For future convenience, we report here the RPA equations for the polarization 
propagator U(k,cu), which is a function of the momentum q and energy uj: 

Tl(q,uj) = U (q,uj) + v c (q)U (q,uj)U(q,uj) . (1) 

Here v c (q) = 4ire 2 /q 2 is the proton-proton Coulomb interaction and Il (g, uj) is 
the free polarization propagator: 

MkM = ^G (|k + qU W))G (kV) (2) 

and G (k, uj) is the free single particle Green' s function 

K ' E-E k + ir] E-E k -irj' K ' 

where E^ is the single particle energy, Uf the Fermi momentum and 9(x) the 
Heaviside step function, which equals 1 for x > and is zero otherwise. The 
excitation spectrum is thus determined by the equation Re(l — v c (k)Uo(k, uj)) = 
0. The explicit expression for the free polarization in the non-relativistic limit 
can be found in textbooks and is usually referred to as the Lindhard function. 
The upper branch appearing in Fig. 2 is called the plasmon excitation, which 
is characterized by a finite energy at vanishing momentum and has a simple 
classical interpretation [16]. The lower branch is actually strongly damped and 
does not correspond to a real excitation which is able to propagate. This point 
will be further elaborated in the sequel. From the Figure it is clear that in any 
case there is a maximum value of the momentum above which the spectrum 
stops and no excitation is possible. This is mainly a quantal effect, which is not 
present in the semi-classical approximation of RPA, which is nothing else that the 
Vlasov approximation. This is illustrated in the same Fig. 2, where the RPA and 
the Vlasov approximation (dashed lines) spectra are compared. Apart from the 
momentum cutoff, one can see that the spectra are quite similar. The results for 



the relativistic case are analogous [15], except that the "thumb" becomes thinner 
and thinner as the Fermi momentum becomes more and more relativistic [17]. 
In the ultra-relativistic case the two branches become extremely close near the 
cutoff, and the difference between RPA and Vlasov spectrum appears smaller and 
smaller, apart again for the presence of the momentum cutoff. In the following, 
all the calculations will be performed in the Vlasov approximation for the three 
fluids, protons, electrons and muons. 

Having summarized the well established results for the excitation spectrum 
with the pure Coulomb interaction within the jelly model, let us consider the 
system of protons, electrons and muons, treating the three components on equal 
footing. The RPA equations for this case are simply a three by three system, 
where the coupling is again provided by the Coulomb interaction between protons, 
muons and electrons: 



The polarization propagator rp fe is now a three by three matrix. In Eq. (4) the 
index i, j, k stands for electrons (e), muons (/i) and protons (p). These equations 
describe the motion of a three-component fluid. It has to be stressed that no 
exchange is included, i.e. only density- density correlations are considered. Since 
both electrons and muons are expected to be faster than protons, they are able 
to follow their motion and are more effective in screening the proton-proton in- 
teractions. The total spectrum is now formed by several branches, but only one 
is expected to have the characteristic property of a finite excitation energy at 
vanishing momentum, and this will be identified with the plasmon mode of the 
proton-muon-electron coupled system. The original proton plasmon mode be- 
comes now a sound-like mode, since, due to the electron and muon screening, the 
effective proton-proton interaction becomes of finite range (screened Coulomb 
interaction). This is a well known phenomenon also in classical plasma [16]. 
The screening can be easily seen in the present case by solving Eqs. (4) for the 
proton-proton polarization propagator. After a short algebra one finds: 



Since the unperturbed electron and muon plasmon modes are substantially higher 
than the corresponding proton one, in the small q limit one can expand the free 
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electron polarization propagator at zero frequency and one gets 

(l " n o P ^) WP = U o P , (6) 

which is an RPA equation with a screened Coulomb interaction and therefore of fi- 
nite range. For q < q c the collective mode can be only a sound mode with a typical 
linear dependence on momentum. The inverse A = l/g c is the screening length. 
One finds that q\ = 3[(ujp/v F ) 2 + (u^/v F ) 2 ], where and tup 1 are the electron and 
muon gas plasmon frequencies, respectively, and v F , v F the corresponding Fermi 
velocities. In neutron star matter conditions q c « k F and the screening length 
is much larger than the average distance between protons. Therefore, at increas- 
ing momenta the proton excitation should merge in the proton plasmon mode 
with no screening. Besides that, the static limit we have performed in Eq. (6) 
is only approximately valid, especially at the higher density, where the proton 
plasma frequency is a substantial fraction of the electron plasma frequency, and 
it is more appropriate to solve explicitly Eq. (4). The actual calculations confirm 
these expectations, as illustrated in Fig. 3 for three proton densities correspond- 
ing to total baryon densities po, 2p and 3po- The reported branches are obtained 
by searching for the energies for which the determinant A of the real part of the 
three by three matrix in Eq. (4) vanishes at a given momentum. In the present 
particular case it has a simple expression: 

a = i - „ c (n« + rc + K) (7) 

One observes that all branches go to zero at vanishing momentum, except the 
upper one. The latter is the plasmon mode of the three components system, while 
the " thumb like" shape of the original proton branch in the uncoupled system can 
be identified with the lowest "loop" (color on line). Additional branches appear 
in between these two. For p = p one additional branch approaches closely 
the plasmon one at increasing momentum. It can be considered as the original 
"sound like" electron mode, according to the above discussion for the proton gas 
in the jelly model. It turns out, as we will show below, that it is still mainly an 
electron excitation even in the coupled case. The last additional excitation, which 
looks also a "sound like" mode is reminiscent of a pure muon excitation, which 
however is now mixed with a small electron component. This coupling become 
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stronger at increasing density and the two branches merge in a single one with 
a typical "bending back", see panels for p = 2p and p = 3po, which actually 
marks the smooth transition form almost a pure electron excitation to a mixed 
electron-muon excitation. 

To support and illustrate this type of analysis we have reported in the same 
figure the branches (dashed lines) one obtains when one considers only protons 
and electrons. In this case, we keep the original proton density and substitue the 
electron-muon liquid by a pure electron liquid with p e = p p . From the figure one 
can trace the evolution of the branches when introducing the muon component. 
It has to be noticed that all branches, except the plasmon one, are actually 
damped. The polarization propagators n which appear in Eqs. (4) are indeed 
in general complex. The best way to estimate the damping rate of each branch 
is to calculate the spectral function — Im(IT M (g, a;)), which, in the present case, 
has a proton (i = p), an electron (i = e) and a muon (i=/i) component. This 
quantity is also the basic one for calculating the rate of many physical processes. 
At a given momentum q, as a function of uj this quantity is expected to display 
peaks corresponding to the different branches, whose width gives a measure of 
their damping rate, i.e. each excitation mode is actually a resonance if the width 
is not too large. The plasmon, up to a certain momentum, has no width, and the 
strength function presents a delta function singularity at the plasmon frequency. 
The electron, muon and and proton strength functions are reported in Fig. 4 
at twice the saturation density and at selected values of the momentum q. For 
guidance, a dashed vertical line is reported at the position of each branch. To 
put in evidence the plasmon branch we have reported just a sharp bar at the 
position of its energy. One can see that that peaks are present at the energy 
corresponding to the branches of Fig. 3, even if some slight shifts occur. Such 
shifts are expected, since the presence of an imaginary part has also the effect to 
produce some displacements. 

The relatively narrow peak at lower energy corresponds to the original proton 
plasmon in the uncoupled system. As discussed above, this branch in the coupled 
system is a sound-like mode due to electron and muon screening. For this branch 
the proton component is larger, but electrons and muons follow the proton os- 
cillation, which allows a quite effective screening. However, as the momentum 
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increases the screening becomes less effective, in agreement with the discussion 
above. Besides that, the coupling with electrons and muons produces a substan- 
tial damping of this branch. At higher energy a small broad peak is obtained, 
which involves a substantial fraction of both muons and electrons, but not pro- 
tons. This branch is a coupled oscillations of electrons and muons, which behave 
approximately as a single fluid. It is characteristic of the coupling between elec- 
trons and muons, and would not be present in a system with only protons and 
electrons. In any case, as one can see, the strength associated with this structure 
is quite low. It has to be noticed that it is present also at q — 10 MeV, where in 
principle no branch is present at this energy. The fact that this structure does 
not correspond to the vanishing of the real part of the RPA determinant is an in- 
dication that the corresponding oscillation is over-damped, and indeed it is quite 
broad at this momentum. 

The very broad structure which appears just below the sharp plasmon peak 
is only due to the electrons, and it corresponds to the over-damped sound-like 
branch which is present in a single fluid, below the plasmon, as discussed for 
the proton case in relation with Fig. 2. Its strength is increasing at increasing 
momentum. 

To illustrate and support these considerations and interpretations of the re- 
sults, we compare in Fig. 5 at different densities and for q — 10 MeV the strength 
functions of the proton-muon-electron system with the ones of the two-component 
proton-electron system, where the electron density equals the electron-muon den- 
sity of the three-component system. One can see that the electrons and muons 
move together in screening the proton plasmon, with a total strength not so far 
from the electron strength in the two-component system. As anticipated above, 
the broad electron-muon structure is not present in the two-component system, 
while the broad electron structure below the plasmon peak is mainly unaffected. 
However, at the lowest density one notices a muon peak just above the proton 
sound-like mode and a very weak involvement of the muons in the screening of 
the proton motion. The muon peak corresponds to the branch at "intermediate 
energy" in Fig. 3. This means that at this density the muons are not able to 
follow the proton motion and produce an independent mode, which however has 
not a plasmon-like character since its energy is mainly linear in momentum. 
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Let us analyze in more detail the plasmon mode, which is the highest in energy 
and appears in the spectral function as a delta-function singularity. In fact at the 
plasmon energy all the free propagators 11" are real. Solving the RPA equations, 
one gets: 

uii = ny[i-(iff + iff> c ] 

A + ir] ' 1 ' 

where % ^ j ^ k and rj is the usual infinitesimal quantity demanded by the 

general structure of the response function. The plasmon contribution to the 
spectral function is then: 

-Im(IP) = v c {W*)H{uj - u p )/\A(w p )'\ (9) 

where oj p is the plasmon energy, solution of the equation A (a;) = 0. The ratios 
of the different strengths are then given by 



Im(IP) (n 



ii\2 
) 



(10) 



Inn II") (ui 3 ) 2 ' 

where the expression of Eq. (7) for the determinant has been used. The plasmon 
pole corresponds also to the zero eigenvalue of the RPA matrix which in Eq. (4) 
defines the response function II (a;), i.e. 

{Sij - n"(u p )v?(q))5 Pj = , (11) 

where the eigenvector Sp = (5p e , 5p^, 5p p ) describes the electron, muon and proton 
density fluctuations of the modes. One easily finds 

a 



5pi_ = zz IIq 



Zizs^h > ( 12 ) 



where Zi — ±1 is the charge sign of the ^-component. This result is in agree- 
ment with Eq. (10), since the spectral function gives the square of the density 
fluctuation. 

In Fig. 6 are reported some ratios of the density fluctuations, according to 
Eq. (12), which describes the plasmon mode at different momenta. The main 
characteristic to be noticed is the dominance of the electron component, with 
the muons becoming of some relevance only at the higher densities. The proton 
component is marginal, since protons are unable to follow the rapid motion of 
the negatively charged components. The momentum dependence is quite weak. 
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In the analysis we have carried out, the nuclear interaction and the presence 
of neutrons have been completely neglected. However, the plasmon mode is at 
relatively large frequency, while the nuclear sound (or zero-sound) mode starts 
linearly with momentum. The coupling between the plasmon and the other modes 
is expected to be quite weak, and the properties of the plasmon excitation are 
likely to be essentially unaffected. Furthermore, the plasmon mode is known not 
to be influenced by the possible presence of a superfluid phase [18] . 

The nuclear interaction can affect the lower energy modes, whose strength will 
be then redistributed. In particular the sound-like mode involving the protons 
will be coupled with neutron modes, since it falls in the low frequency part of 
the spectrum, due to the screening effect discussed above. This show the rele- 
vance of the Coulomb interaction for the overall structure of the spectrum. A 
complete analysis of the excitation strength distribution, including both Coulomb 
and nuclear interaction, will be given elsewhere. 

III. CONCLUSION AND PROSPECTS. 

We have analyzed, within the RPA scheme, the plasmon excitations which can 
occur in the homogeneous matter present in the outer core of Neutron Stars. The 
proton, electron and muon components have been considered. Since both electron 
and muon Fermi velocity are larger than the proton one, the Coulomb proton- 
proton interaction is effectively screened and no genuine proton plasmon can exist, 
but instead a sound-like mode is present. The analysis of the spectral functions 
shows that the high frequency plasmon mode is mainly an electron excitation, 
with possibly a substantial muon component at higher density. Except for the 
genuine plasmon mode, all other excitation modes are damped, and in some 
case over-damped. The inclusion of the nuclear interaction and of the neutron 
component will be reported elsewhere. 
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Ph [fm -3 ] 


Y P [%] 


Y e [%} 


Yp [%] 


w 0e [MeV] 


0.16 


3.70 


3.63 


0.07 


6.15 


0.32 


8.57 


5.95 


2.62 


10.25 


0.48 


13.03 


7.98 


5.05 


13.48 



TABLE I: Particle fractions for the baryonic densities under study. The electron plas- 
mon frequency for p e = p p , used as a unit for the plots, is also given. 

Figure captions. 

Fig. 1 - The proton fraction (thick lines) and muon fraction (thin lines) as a 
function of total baryon density in neutron star matter. Microscopic calculations 
(taken as a reference) are compared with results from modern Skyrme forces. 

Fig. 2 - The proton excitation spectra in the jelly model. The energy u>o e is 
the electron plasmon energy at zero momentum, for p e = p p . The full line is 
obtained in the RPA scheme, the long-dashed line corresponds to the Vlasov 
approximation. The dashed-dotted line shows the effect of electron screening on 
the (RPA) proton spectrum. 

Fig. 3 - Dispersion relation in the pep, system (full lines), compared with the pe 
system (dashed lines), for three baryonic densities. The corresponding particle 
fractions are given in Tab. I. The dotted lines correspond to Ui = vpiq. The 
energy Uo e is the electron plasmon energy for p e = p p . 

Fig. 4 - Diagonal elements of the spectral-function matrix {pep: system) , for three 
momenta and pb = 0.32 fm -3 . The vertical lines indicate the energy at which the 
real part of the determinant of Eq. 7 vanishes. For the plasmon the dashed and 
full lines coincide. The energy u 0e is the electron plasmon energy for p e = p p . 

Fig. 5 - Diagonal elements of the spectral-function matrix (pep, system) , for three 
different densities and q = 10 MeV. The energy u 0e is the electron plasmon energy 
for p e = p p . The vertical lines indicate the energy at which the real part of the 
determinant of Eq. 7 vanishes. For the plasmon the dashed and full lines coincide. 
Left: pep, system; right: pe system. 

Fig. 6 - Direction of the eigenvector associated with the plasmon mode. 
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